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1. FORMAT 

Text should be produced within the dimensions 
shown on these pages: each column 7.5 cm wide 
with 1 cm middle margin, total width of 16 cm 
and a maximum length of 20.2 cm on first pages 
and 21 cm on second and following pages. The 
I^T^^X document style uses the maximal stipu- 
lated length apart from the following two excep- 
tions (i) M^TpC does not begin a new section di- 
rectly at the bottom of a page, but transfers the 
heading to the top of the next page; (ii) I^T[t;X 
never (well, hardly ever) exceeds the length of the 
text area in order to complete a section of text or 
a paragraph. 

1.1. Spacing 

We normally recommend the use of 1.0 (sin- 
gle) line spacing. However, when typing com- 
plicated mathematical text BT[t;X automatically 
increases the space between text lines in order 
to prevent sub- and superscript fonts overlapping 
one another and making your printed matter il- 
legible. 

1.2. Fonts 

These instructions have been produced using 
a 10 point Computer Modern Roman. Other 
recommended fonts are 10 point Times Roman, 



* Footnotes should appear on the first page only to indi- 
cate your present address (if different from your normal 
address), research grant, sponsoring agency, etc. These 
are obtained with the Xthanks command. 



New Century Schoolbook, Bookman Light and 
Palatino. 



2. PRINTOUT 

The most suitable printer is a laser printer. A 
dot matrix printer should only be used if it pos- 
sesses an 18 or 24 pin printhead ("letter-quality"). 

The printout submitted should be an original; 
a photocopy is not acceptable. Please make use 
of good quality plain white A4 (or US Letter) 
paper size. The dimensions shown here should be 
strictly adhered to: do not make changes to these 
dimensions, which are determined by the docu- 
ment style. The document style leaves at least 
3 cm at the top of the page before the head, which 
contains the page number. 

Printers sometimes produce text which con- 
tains light and dark streaks, or has considerable 
lighting variation either between left-hand and 
right-hand margins or between text heads and 
bottoms. To achieve optimal reproduction qual- 
ity, the contrast of text lettering must be uniform, 
sharp and dark over the whole page and through- 
out the article. 

If corrections are made to the text, print com- 
pletely new replacement pages. The contrast on 
these pages should be consistent with the rest 
of the paper as should text dimensions and font 
sizes. 
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Table 1 



Biologically treated effluents (mg/1) 




Pilot plant 




Full scale plant 






Influent 


Effluent 


Influent 


Effluent 


Total cyanide 


6.5 


0.35 


2.0 


0.30 


Method-C cyanide 


4.1 


0.05 




0.02 


Thiocyanide 


60.0 


1.0 


50.0 


< 0.10 


Ammonia 


6.0 


0.50 




0.10 


Copper 


1.0 


0.04 


1.0 


0.05 


Suspended solids 








< 10.0 
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Ritcey, Tailings 
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Amsterdam, 
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3. TABLES AND ILLUSTRATIONS 

Tables should be made with WT^: illustra- 
tions should be originals or sharp prints. They 
should be arranged throughout the text and 
preferably be included on the same page as they 
are first discussed. They should have a self- 
contained caption and be positioned in flush- 
left alignment with the text margin within the 
column. If they do not fit into one column 
they may be placed across both columns (us- 
ing \begin{table*} or \begin{f igure*} so that 
they appear at the top of a page). 

3.1. Tables 

Tables should be presented in the form shown 
in Table |l|. Their layout should be consistent 
throughout. 

Horizontal lines should be placed above and be- 
low table headings, above the subheadings and at 
the end of the table above any notes. Vertical 
lines should be avoided. 

If a table is too long to fit onto one page, 
the table number and headings should be re- 
peated above the continuation of the table. For 
this you have to reset the table counter with 
\addtocounter{table}{-l}. Alternatively, the 
table can be turned by 90° ('landscape mode') 
and spread over two consecutive pages (first 
an even-numbered, then an odd-numbered one) 
created by means of \begin{table} [h] without 
a caption. To do this, you prepare the table as 
a separate I^Ti?;X document and attach the tables 
to the empty pages with a few spots of suitable 
glue. 



3.2. Line drawings 

Line drawings should be drawn in India ink 
on tracing paper with the aid of a stencil or 
should be glossy prints of the same; computer pre- 
pared drawings are also acceptable. They should 
be attached to your manuscript page, correctly 
aligned, using suitable glue and not transparent 
tape. When placing a figure at the top of a page, 
the top of the figure should be at the same level 
as the bottom of the first text line. 

All notations and lettering should be no less 
than 2 mm high. The use of heavy black, bold 
lettering should be avoided as this will look un- 
pleasantly dark when printed. 

3.3. Black and v^rhite photographs 

Photographs must always be sharp originals 
(not screened versions) and rich in contrast. 
They will undergo the same reduction as the text 
and should be pasted on your page in the same 
way as line drawings. 

3.4. Colour photographs 

Sharp originals (not transparencies or slides) 
should be submitted close to the size expected 
in publication. Charges for the processing and 
printing of colour will be passed on to the 
author(s) of the paper. As costs involved are per 
page, care should be taken in the selection of size 
and shape so that two or more illustrations may 
be fitted together on one page. Please contact 
the Technical Editor in the Camera-Ready Pub- 
lications Department at Elsevier for a price quo- 
tation and layout instructions before producing 
your paper in its final form. 
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Figure 1. Good sharp prints should be used and 
not (distorted) photocopies. 



in parentheses flush with the right-hand margin 
of the text and level with the last line of the equa- 
tion. For multi-line equations, use the eqnarray 
environment. For complex mathematics, use the 
Am.S-WT^ package. 
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References should be collected at the end of 
your paper. Do not begin them on a new page 
unless this is absolutely necessary. They should 
be prepared according to the sequential numeric 
system making sure that all material mentioned 
is generally available to the reader. Use \cite 
to refer to the entries in the bibliography so that 
your accumulated list corresponds to the citations 
made in the text body. 

Above we have listed some references according 
to the sequential numeric system ||l|-^ . 



Figure 2. Remember to keep details clear and 
large enough. 



4. EQUATIONS 

Equations should be flush-left with the text 
margin; I^TJ?;X ensures that the equation is pre- 
ceded and followed by one line of white space. 
M'eX provides the document-style option f leqn 
to get the flush- left effect. 

H^piiu) = E^^\co)S^p + {a\W^\l3) (1) 

You need not put in equation numbers, since 
this is taken care of automatically. The equation 
numbers are always consecutive and are printed 
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Modelling the Quark Determinant in Full QCD Simulations 

A. Duncan'^jE. Eichten^,and H. Thacker'^ 

'^Dept. of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260 
^ Fermilab, PO Box 500, Batavia, IL60510 

■^Dept. of Physics, University of Virginia, Charlottesville, VA 22901 

The computational requirements and dynamics of Monte Carlo simulations of unquenched QCD incorporating 
the infrared quark eigenmodes (up to ~ ^qcd) exactly and UV modes via a loop representation are discussed. 
The accuracy of such a loop representation is studied for a variety of lattice volumes and quark masses. The 
method has been successfully applied for lattices up to 10^x20 at a ~0.17F with improved (clover) action, and 
allows simulations at or near kappa critical. 
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1. Splitting the Quark Determinant 

The essence of the truncated determinant ap- 
proach ||l| to unquenched QCD lies in the real- 
ization that the infrared part of the quark deter- 
minant (specifically, the determinant ^{A) of the 
hermitian operator H = 75(1/ (A) — m)) can be 
gauge-invariantly split off, leaving an ultraviolet 
part which is accurately fit by a linear combi- 
nation of a small number of Wilson loops. The 
eigenvalues of H 

(1) measure quark off-shellness (for A = 0, Xi 

(2) are gauge-invariant, Xi{A) ~ Xi{A^). 

Thus we can write V{A) = Vnz{A)Vuv{A), 
where the infrared part Vxui^) defined as the 
product of the lowest Nx positive and negative 
eigenvalues of H, with |Ai| < Aqcd (typically, — 
300-400 MeV). This cutoff is chosen (a) to include 
as much as possible of the important low-energy 
chiral physics of the unquenched theory while (b) 
leaving the fluctuations of In 'Dxtz of order unity 
after each sweep updating all links with the pure 
gauge action. This ensures that the acceptance 
rate is sufficiently high when the infrared deter- 
minant only is used in the accept/reject stage of 
the procedure. The crucial point is that it is pos- 
sible to achieve both (a) and (b) on fairly large 
lattices (up to physical volume ~ 20 F**) as well 
as at kappa values arbitrarily close to kappa crit- 
ical. 



2. Efficient Computation of the IR spec- 
trum in QCD4 

The following Lanczos procedure allows us to 
extract the needed infrared eigenvalues of H rel- 
atively rapidly: 

(1) Starting from an initial vector vi , an orthonor- 
mal sequence vi, V2, V3, ..vnl is generated by the 
standard recursion: 

with the constants at, (3k determined from over- 
laps of generated vectors. In the basis of the Vi, H 
is tridiagonal. The corresponding real symmetric 
tridiagonal matrix T^^ has the au on the diago- 
nal and the (3k on the sub (and super) diagonal. 

(2) A CuUum-Willoughby sieve Q is used to iden- 
tify and remove spurious eigenvalues. 

(3) The remaining "good" eigenvalues converge 
most rapidly in the least dense part of the spec- 
trum, in particular, in the needed infrared por- 
tion. This remains true even at kappa critical ( in 
marked contrast to the matrix inversions needed 
in HMC simulations, for example), where the den- 
sity of eigenvalues near zero is still small. The sta- 
bility and accuracy of the converged eigenvalues 
has been checked extensively by gauge transform- 
ing the gauge field. 

(4) The diagonalization of the T^^ matrix (typ- 
ically, of order 10,000 in the QCD case) can 
be completely parallelized using the Sturm se- 
quence property H of tridiagonal matrices in 
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which nonoverlapping parts of the spectrum are 
independently extracted by a bisection procedure. 

3. Fitting the UV modes- Loop represen- 
tations of the Quark determinant 

The effective gauge action generated by inter- 
nal quark loops is a gauge-invariant functional, 
which can be evaluated explicitly in a hopping 
parameter expansion: 

\nV{A)/V = 288k'* ^ Li -f 2304k^ ^ L2 

+ 4608k^^L3 + 1536k6^L4 + .. 

where V is the lattice volume, Li a generic pla- 
quette, and ^2,3,4 are 6 link loops with link direc- 
tions {i,j,k,~j,-i,-k), and 
{i,j,k,-i,~j,~k) (1 < i,j,k < 4). Unfortu- 
nately, this hopping parameter expansion is use- 
less Q except for extremely heavy quarks (k — > 
0), due to the contribution of large loops. In- 
stead, large loops may be cut off by removing 
the IR modes |Ai| < Aqcd- Now the expansion 
converges much more quickly, and we may write 
InVuvi^) = yJ2^ii^) where the sum involves 
only a small number of loops (as we shall see, typi- 
cally less than 10) and the coefficients Ci are deter- 
mined nonperturbatively. Of course, to compute 
T^uv , we need the complete spectra for an ensem- 
ble of configurations. The computational cost for 
extracting a complete Dirac spectrum via Lanc- 
zos is large but manageable, as this calculation 
need only be done for a limited number of decor- 
related configurations. For example, the complete 
Dirac spectrum for a 10^x20 lattice has 240,000 
eigenvalues and requires about 500,000 Lanczos 
sweeps, equivalent to about 1 400-Mhz-Pentium- 
week. The final spectrum can be checked with 
analytic spectral sum rules which give the sum of 
powers of the eigenvalues as explicit functionals 
of small loops. 

The results of a fit of In Vjjv to a linear combi- 
nation of Wilson loops on an ensemble of 10'^x20 
lattices at beta=5.7 (with clover improvement) 
and K=0.1425 is shown in Fig.l The configura- 
tions were generated including the truncated de- 
terminant InP/fl and therefore already contain 
the exact low-energy chiral physics. The fit is 
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Figure 1. Fit of UV Determinant to Loop Action, 
K=0.1425 



very good once 4, 6 and 3 figure-8 8-link opera- 
tors are included in the fit . These results sug- 
gest that the full determinant can be accurately 
modelled by computing the low eigenvalues ex- 
actly and including the remaining high modes 
via an approximate loop action. In any event, 
we expect that the UV fluctuations affect pri- 
marily the scale of the theory, while for the low 
energy spectrum, quark off-shellness is limited 
to about Aqcd and dimensionless mass ratios 
should therefore be largely insensitive to Vjjv- 

4. MonteCarlo Dynamics for QCD4 Simu- 
lations with the Truncated Determinant 

In the truncated determinant approach to 
QCD4 Ql, unquenched configurations are gener- 
ated by the following algorithm: 

1. Update the gauge configuration with the 
pure gauge action, using a procedure com- 
pliant with detailed balance. We have used 
Metropolis link updates applied to a ran- 
domly chosen block of noninterfering links 
to ensure detailed balance while maintain- 
ing parallelizability of the computation. 

2. Apply a metropolis accept/reject criterion 
based on the the effective quark action 

= NplndeWiRiA) {Np = 2) (1) 
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Figure 2. Determinant relaxation, k;=0.1440 



Typically, we use an IR cutoff for the truncated 
determinant corresponding to a gauge-invariant 
eigenvalue of 300-400 MeV. This procedure leads 
to tolerable acceptance rates on lattices of fairly 
large physical volume - we have explored systems 
up to 20 -F*^, while controllable finite size errors 
in electromagnetic fine structure lattice studies 
with a long range massless U(l) field require lat- 
tice volumes >6 F^. The simulations in progress 
are on three different lattice sizes: 

(1) 10'^x20 lattices at /3=5.7 (clover improved) 
at K=0.1415, 0.1425, 0.1436 and 0.1440 (the 
last value being very close to kappa critical). 
This corresponds to a physical volume of roughly 
(1.7F)3x3.4F. 

(2) 6'' lattices at an effective (3 of 4.5, but us- 
ing the O(a^) improved gauge action of Afford et 
al, Q , at kappa critical. 

(3) 8^ lattices at an effective (3—4.5 (O(a^) im- 
proved), again at kappa critical. 

For the 10^x20 lattices we have checked that 
the truncated determinant simulations succeed 
in equilibrating the configurations, and that 
they decorrelate reasonably rapidly subsequent 
to equilibration. The equilibration can be stud- 
ied by looking at the relaxation of S'quark from 
the quenched value corresponding to the starting 
configuration. Even a,t k = Kc, the 10^x20 lat- 
tices equilibrate after a few hundred sweeps (see 
Fig. 2). 

To measure decorrelation we have calculated 
the autocorrelation of the pion propagator at vari- 
ous time separations, keeping configurations sepa- 
rated by 20 steps of the basic algorithm. A typical 




config # 

Figure 3. Decorrelation of pion propagator, 
K=0.1425 



example, from the 10^x20 runs with k=0.1425, is 
shown in Fig. 3, where the pion correlator at time 
slice 6 is seen to be effectively decorrelated after 
about 30 steps of the algorithm. 

Cases (2,3) represent physically large lattices 
(33 F"* and 105 F'* resp.) at kappa critical and 
display critical slowing down (several thousand 
sweeps are needed to equilibrate the configura- 
tions). However the lattice sizes are small: new 
configurations can be generated and the trun- 
cated determinant computed quite rapidly (~ 20 
minutes for the 6"^ lattices, two hours for the 
8^ lattices, on a Pentium 400 Mhz processor). 
These lattices will allow a detailed study of string- 
breaking and other unquenched dynamical effects 
(see talk of Eichten, this conference). 
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